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Abstract 

We provide a geometrical interpretation for the series of transformations used by Sakovich to 
map the third-order nonlinear evolution equation obtained by Chou and Qu to the mKdV equation. 
We also discuss its bi-Hamiltonian integrability as well as integrable equations associated with this 
system. 
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1 Introduction: 



Starting with the work of Hasimoto pQ the study of the connection between integrable models 
and the dynamics of space curves has drawn a lot of interest. Motion of these curves and the 
associated integrable models were considered in Euclidean space E 3 [2] [3j 13] , S 2 and S 3 spaces [5] 
and Minkowski space [6j [7]. However, integrable equations associated with the motion of plane 
curves are more scarce [SJ E]. This has motivated Chou and Qu [TUJ [TT] to study motions of 
plane curves in other background geometries, replacing the Euclidean geometries by Klein ones. 
As a result they have shown that KdV, Harry-Dym and Sawada-Kotera hierachies, among other 
nonlinear equations, arise in an natural way from the motions of plane curves in affine, centro- 
affine and similarity geometries. In this process they also found a new equation associated with 
the mKdV equation, namely (see [ID], last equation on page 31), 

ut = -(iu xx +u)- 2 ) x . (1) 

Chou and Qu did not investigated if this equation (from now on the CQ equation) arises from an 
AKNS- or the WKI-scheme of inverse scattering transformation, its integrability was investigated 
by Sakovich [T2] . Through a chain of Miura-type transformations he related equation ([!]) with the 
mKdV equation. Firstly, (P) can be transformed to 

v t = l y2 (( v2 )^ + v2 ) x > ( 2 ) 

by the transformation 

(x,t,u(x,t)) i->- (x,t,v(x,t)) : v = -(u xx + uy 1 . (3) 
Then we perform two successive transformations 

(y,t,w(y,t)) (->■ (x,t,v(x,t)) : x = w , v = w y , 

(y,t,w(y,t)) ^ (y,t,z(y,t)) : z = w y , (4) 
to obtain the mKdV equation 

3 2 / r- \ 

z t = z yyy + ~ z z y , (5) 

where z(y,t) = v(x,t). Also, in [12J a zero curvature representation with an essential parameter 
was obtained as well as the following second-order recursion operator 

r = — - — d — - — {d + d- 1 ). (6) 

u xx + u u xx + u 

It is well known that integrable equations also possess a bi-Hamiltonian structure which yields 
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a recursion operator in the form T>-iDy . So, going backwards Ren and Alatancang [13] proved the 
hereditary property of © and by a decomposition of it obtained the bi-Hamiltonian structure as 
well as a hierarchy of equations, associated with ([1]). 

In this paper we point out, in Section 2, that the series of transformations ([3]) and dH) mapping 
the CQ equation ([1]) to the mKdV equation © has a nice geometrical interpretation which can 
be used to justify connections between other equations. In Section 3 the bi-Hamiltonian structure, 
and consequently the recursion operator R, is derived directly from a Lagrangian representation 
of (pQ). In fact, the original derivation of Ren and Alatancang for this bi-Hamiltonian structure 
relies on the knowledge of the bi-Hamiltonian structure of a system related to ([T]) by a Miura map 
(obtained previously by Olver and Rosenau |14| ) and of the recursion operator ©. Therefore, 
their derivation is more a map between structures than a derivation from first principles. For 
completeness, in Section 4, we obtain the hierarchy of equations, some symmetries and equations 
related to ([T]). We discuss the relation among these equations. The conclusions are presented in 
Section 5. 



2 Plane Curve Motion: 



Let be a closed smooth curve in the plane, parametrized by an arbitrary parameter a, where 
a G [0,1], and let r(a,i) represent the position of a point on the curve at the time t. Along the 
curve the arc-length is s(a,t) = y g(a' ,t) da' and we can use s as a parameter as well (this is 
called Lagrangian description of the curve). The metric on the curve is g = r a ■ r a and we have 
that d/da = ^fgdjds. At a point on the curve the unit tangent and normal vectors, defined by 
t(s,t) = r s and n(s,i) = ^r ss , respectively, satisfy the Serret-Frenet equations 




where k(s, t) is the curvature. The dynamics of the points of the curve are specified by 

r t = Fn + Gt , (8) 

where the normal and tangential velocities F and G are functions of the curvature k. Now we make 
the further assumption that the perimeter L = j> ds of the closed curve remains constant in time, 
i.e., we assume an isoperimetric plane curve motion. Therefore, from the metric evolution 

g t = 2r a ■ r ta = 2g (G s - kF) , 



we get 



— = J y /g(a , ,t)da' = j (G s - kF) ds , 
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and dL/dt = if G s = kF. This implies that the arc- length and time derivatives commute, 
[d s ,dt] = 0. The time evolution of the tangent and normal vectors gives 



F s + kG 

-F. -kG 



ia 2 (F s + kG) 



n 



and the compatibility condition between ([7|) and ([9]), 



ts \ / st 

yields the following nonlinear evolution equation for the curvature 

Kt = F ss + Ks d; 1 ( K F) + K 2 F , 
which can also be written as Kt = RF with 

R = dt + KgdT 1 !^ + k 2 . 



(9) 



(10) 
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Figure 1 

If we choose F = —k s in (|10p we obtain the mKdV equation, 

_ 3 2 

Kt- Ksss 2 nn s , 



(12) 



and we recognize (jlip as the mKdV recursion operator. However, at this point let us change the 
curve description and parametrization. For a curve with positive curvature (uniformly convex) we 
introduce the support function h(0,t), defined as 



h(9,t) 



r • n 



(13) 
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where 9, called the normal angle of the curve, is the angle between the unit normal n and the positive 
a>axis [15]. In the following let r = r(9,t) and from the Figure 1 we have n = — (cos#,sin#) and 
t = (— sin 6, cos 9). For the support function h(9,t) we obtain 

h = x cos 9 + y sin 9 
hg = —x sin 9 + y cos 9 , 

and r = (x, y) can be obtained form h using 

x = h cos 9 — hg sin 9 
y = h sin 9 + hg cos 9 . 

From this last equation we obtain 

M = hgg + h , 

and since d/ds = \rg\~ 1 d/d9 the Serrat-Frenet formula 

dt i dt 

ds = d9 = KU 

yields a neat expression for the curvature in terms of the support function 

' (14) 



hgg + h 

Finally, from ([8]) the normal speed of r is given by n • = — k s and therefore 

h t = k s = Irel" 1 /^ = ^(/« 2 )e 

= \ {{hgg + h)- 2 ) e , (15) 

which is the equation (TjQ). An evolution equation for the curvature, in terms of the normal angle 9, 
instead of the arc-length as in ()12p . can be obtained. From (|14[) 

K t = -K 2 (hget + h t ) , 

and after using (fT5|) 

«t = -5« 2 ((« 2 )w + « 2 ) fl - (16) 

It is well know that a curve is determined by its curvature up to a rigid motion [T2], therefore, 
there is a formal equivalence between (|16p and f)15|) (and between (|8|) and (|12p ) and the connection 
is made by (|14p . In other words the equation (|15p and the curvature integrable equation (|16p are 
equivalent. Now we are in position to justify geometrically the transformations Q and if we 
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make the identifications h(9,t) o u(x,t), n(s,t) o z(y,t), K(9,t) -R- v(x,t), 9 x and soy 
and recognize the equivalence between the equations (fT5|) o ([I]), (fT2|) © and (fT6j) <-> ([2]). The 
transformation ([3]) is simply (|14p and the transformation (J3|) is the inverse transformation from 
(|12|) to (|16|) (projection (|13|) and change of parametrization from s to 9). The recursion operator 
([6]) defines the curvature evolution under arc-length conserving dynamics and can be obtained from 
(fTT|) using the identifications given above. 



3 Bi-Hamiltonian Structure: 

Let us study the integrability of (|15|) from a Hamiltonian point of view. Since the curvature plays 
a prominent role let us write (fT5j) as 

ht = \(^)e , (17) 
and where k satisfies (|16p which can also be rewritten as 

(K- l ) t = l -[(K 2 ) ee + K 1 ) e . (18) 

Let us observe that the basic field is h and that k is a placeholder used to make expressions 
more compact (see [E]). The equations (|16p and (|18p are not being viewed as nonlinear evolution 
equations but just as expressions for the time derivative of the placeholder k, however, they can be 
interpreted as the Casimir equation for the so called modified compacton hierarchy (see equations 
(27) and (28) in [2]). Also, we will consider ^-periodic solutions, therefore, 9 £ [0, 2tt] and h(9, t) = 
h(9 + 27T, t). We denote by d = dg the differential operator with respect to 9 with skew-adjoint 
inverse d^ 1 = 87 ■ Whenever the nonlocal operator d~ l is used we have dd~ l = d~ 1 d = I, 
= — d and (<9 -1 )^ = —d~ l . We define the anti-derivative d~ l as d9'. 

Equation (|17p can be obtained from a variational principle, 5 J dtd9C, from the Lagrangian 
density 

£ = K-htd- 1 ^ 1 . (19) 

This is a first order Lagrangian density and we can use the Dirac's theory of constraints [T7] to 
obtain the Hamiltonian and the Hamiltonian operator associated with (|19p . The Lagrangian is 
degenerate and the primary constraint is obtained to be 

<S> = Tr + h e + (d- 1 h) , (20) 

where it = dC/dht is the canonical momentum. The total Hamiltonian is 

H T = j d9 {irh t - C + A$) = J ' d9[-K + \(ir + h e + (d _1 /i))] , (21) 
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where A is a Lagrange multiplier field. From the canonical Poisson bracket relation 

{h(9),n(9')} = 5(9-9'), (22) 
with all others vanishing, we require the primary constraint to be stationary under time evolution, 

me),H T } = o, 

to determine the Lagrange multiplier field A in (|2ip and to find out that the system has no further 
constraints. 

The canonical Poisson bracket relation (|22p yields 

K(9,9') = {$(9),$(9')} = (d e + dg 1 ) 5(6 - 9') - (d e , + de,- 1 ) 5(9' - 6) , (23) 

and we find that the constraint (|2Q|) is second class. The Dirac bracket between the basic variables 
is 



{h(9),h(9')} D = {h(9),h(9')} - J d9 1 d9 2 {h(9), <S>(9 1 )}J(9 1 ,9 2 ){$(9 2 ),h(9')} = J(9, 9') , 
where J is the inverse of the Poisson bracket of the constraint ()23|) . 

J d9" K(9, 9")J(9", 9') = 5(9 - 9') . 

From this last equation we get 

2(8 + d- 1 ) J(9,9') = 5(9-9') , 



or 



J(9,9') =V l 5(9-9') , 

where 
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We now set the constraint (I20p strongly to zero in (121 h to obtain 



V 1 = l -(d + d- 1 Y 1 . (24) 



H 2 = - d9K, (25) 



and the equation f)17|) can be written in the Hamiltonian form as 

t= 1 ~5h' ' 



From (|18p it is straightforward to show that the charge 

Hi = j d9^ (26) 
is also conserved. Therefore, the equation ()17p can be written in the Hamiltonian form as 

* = 2 ~5h ' 

where we have defined 

£> 2 = k8k . (27) 
This Hamiltoninan structure is manifestly skew symmetric. Using the expansion 

oo 

(d + d- 1 )' 1 = Y,(-l) n d~ [2n+1) (28) 

n=0 

it is easy to show that T>\ is also skew symmetric. Jacobi identity for these structures as well as 
their compatibility follows from standard method of prolongation [T5]. We can construct the two 
bivectors associated with the two structures as 



®v 2 = ~ J d9 {eft AV 2 <ft} = ~ J deK 2 <f>A<f>, 



-n- 1 



where (ft is the univector corresponding to the one-form dh. Using the prolongation relations, 

prv Vl<j> (K 2 ) = -t^cfto , 

pr v v ^{k 2 ) = -2k 3 [(8 + d- 1 ) K(K(ft)e] e , (29) 
we show that the prolongation of the bivector <3?x> 2 vanishes, 

prvv^^vz) = , 
implying that T>2 satisfies Jacobi identity. Using ([29jh it also follows that 

pri^($x> 2 ) + priW (^i) = > 

showing that T>\ and T>2 are compatible. Namely, not only are T>i,T>2 genuine Hamiltonian struc- 
tures, any arbitrary linear combination of them is as well. As a result, the dynamical equation (|17p 
is bi-Hamiltonian and, consequently, is integrable [18|. I19J . Also, the bi-Hamiltonian structures (I24j) 
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and (|27p provide us with a natural recursion operator defined by 

R = V 2 V± 1 , 



or 



R = 2 K dn(d + d- 1 ) , 
which is exactly the recursion operator R, given by ([6j), obtained by Sakovich [12 



4 Hierarchy of Equations: 

From the bi-Hamiltonian structure obtained in the previous section we can naturally define a 
hierarchy of commuting flows by 



htn =K n [h]=V l 5 ^=V 2 5 ^, n=l,2,3,.... 



(30) 



For n = lwe obtain (fTT)) . From (j30|) we have 



5H n +\ „±5Hr, 



5h 



(31) 



where 

R\ = (V 2 V^ = 2(<9~ 1 + d)ndK 
is the adjoint of R. Using (131 p recursively we obtain an infinite set of conserved Hamiltonians 

Hi = 



d0- 



Ho 



H 3 



Ha 



d0K 



d6 (k 3 - 4,KK 2 e ) 



-24 J d6 (3k 5 + 20k 3 K 2 ee + 15k 4 kqq + n 4 K deee ) 



-2 I ,10 

28 



5 j 35 a ( 13 a 2T r: o 

16 12 V 120 20 e6 



-5^2 



21 



15 



4 2 
K K/i 



21 4 3 16 6 
2 ee 45 
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(32) 



These charges can also be obtained using the fact that 



H n = TrR— , n = 1,2,3,... , 



where "Tr" is the Adler's trace [20] . We also obtained (|32p in this way using our program PSEUDO 

eh. 



The corresponding first flows associated with (|32j) are 
h ti = \ (« 2 )e ' 

h t2 = ( ^K 2 + 2K 3 K Sg + K 2 K% 



h t3 = ( ^k 6 + 15k 4 Kg + 10k 5 Kee + ISk 4 ^ + 2Ak 4l K0K000 + 4k 5 K^e + 2k 2 k# + 32k 3 k^k^ 



(33) 



Since every symmetry of an integrable model defines another integrable model we can go further 
in the study of the hierarchy of equations of our system through a symmetry study of Eq. (|17j) . 
Using the Lie's algorithm (assisted by the computer algebra system program GeM [22] ) we have 
obtained the following point and first higher order symmetry generators in evolutionary form 



X, 



Xo 



X, 



X, 



X, 



hg 



d_ 

dh 



sm( 



d_ 

dh 



h 

+ 3 



d_ 
dh 



Of course the flow associated with X4 is in the hierarchy (|33p . The flow associated with X2 can be 
obtained if we allow in (|30p the value n = 0. In fact, from 



h 



to 



V 



~Sh ' 
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where 5H\/8h = 1, and after using (|28p and (d 1 1) = 8 n /nl we obtain 



^ = \ (d + d-y 1 • i = i£(-i) n (^ (2n+1) i) = isine . 

n=0 

From the flow Xi we take ht = ho as a seed equation and the recursion procedure /ij = (i? -1 )" - ^ 
generates for n = 1 

h t + = - + /i ee )(/i 2 + h 2 e )] e . 

This is the dual counterpart of the mKdV equation considered by Olver and Rosenau (equation 
(25) in [14]). In this same work we recognize (up to multiplicative constants) equation (|16p as the 
Lagrange transform of the mKdV equation (equation (28) in [E]). Prom f)16|) (or (|18p ) if we make 
the identification k = 1/pwe obtain 

Pt = \(p- 2 )eee + \{p- 2 )e, (34) 

called Casimir equation for the modified compacton hierarchy (up to multiplicative constants it is 
the equation (27) in [14]). It is implicit in p3] (due to the role played by the mKdV equation) that 
(|34"j) and the CQ equation (fT5j) are related by the Miura transformation p = u + u xx . Using this 
transformation the bi-Hamiltonian structure of ()34|) follows from the bi-Hamiltonian structure of 
the CQ equation derived in Section 3 and they yield 

V 1 = d + d 3 , 
T>2 = dpd~ 1 pd , 

The input of this bi-Hamiltonian structure was used in [13] for the derivation of the bi-Hamiltonian 
structure of the CQ equation. In fact they transformed T>±, H2 back to the CQ equation variables 
and the second Hamiltonian structure was obtained by a factorization of the recursion operator 
®. 

5 Conclusion: 

In this paper, the transformations used by Sakovich to map the CQ and mKdV equations into each 
other were shown to have a nice geometrical interpretation. Namely, they follow from a isoperimetric 
curve motion in the Euclidean plane. Also, the bi-Hamiltonian structure of the system was obtained 
without any knowledge of its recursion operator or other known result of equations associated with 
the system under consideration. Finally, we have derived a hierarchy of equations, symmetries and 
equations associated with the CQ equation system which have appeared in the literature previously. 



H- 



dOp- 1 , 

Hl = - l -Jde { P - 3 -Ap- 5 p 2 e ) 
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